extern double *u, J, M, L, dx, alpha,beta;
extern int steps,test_n;
extern void get_u();
